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A linear-quadratic-Gaussian (LQG) compensator design procedure is proposed 
for the DSS-13 antenna . The procedure is based on two properties . It is shown that 
tracking and flexible motion of the antenna are almost independent (the separation 
property) , As a consequence , compensators for the flexible and tracking parts can 
be designed separately. It is shown also that the balanced LQG compensator’s effort 
is evenly divided between the controller and the estimator. This allows a minimiza- 
tion of the compensator order, which is important for implementation purposes . 
An efficient compensator reduction procedure that gives a stable low-order com- 
pensator of satisfactory performance is introduced. This approach is illustrated 
with a detailed compensator design for the DSS-13 antenna. The implementation 
of this compensator design requires an update of the antenna model. 


I. Introduction 

The linear quadratic controller for the DSS-14 antenna 
was designed by Alvarez and Nickerson [1], and a linear 
quadratic controller for the DSS-13 antenna was designed 
by Gawronski [2]. The design method presented in this 
article extends the results obtained in [2] for the case when 
full-state feedback is not available. 

The development of new high-performance controllers 
for the DSN antennas is a current priority. The existing 
proportional and integral (PI) controllers satisfy the re- 
quirements for X-band (8.4-GHz) tracking; they remain 
simple, robust to parameter variations, and do not re- 
quire detailed knowledge of the antenna dynamics. How- 
ever, due to the recent pointing requirements for Ka-band 
(32 GHz), new performance requirements for the antenna 
controllers have emerged. The PI controllers have reached 
their performance limits, therefore a new generation of 
controllers has to be designed and developed. Also, in 


order to improve controller performance, more sophisti- 
cated and accurate antenna models have to be developed. 
As a rule, the better the knowledge of the plant dynam- 
ics, the better the performance that can be achieved by 
the controller. The recently developed antenna models [3] 
are accurate enough to give an opportunity to improve 
tracking performance. The models allow simulation of si- 
multaneous tracking in azimuth and elevation, and include 
antenna flexible deformations up to 10 Hz. 

Among the family of the model-based controllers, the 
linear- quadratic- Gaussian (LQG) compensator has been 
chosen for the DSN antennas because it is commonly 
known to be performing well in industrial applications. 
The LQG compensator consists of a controller and an esti- 
mator. The controller drives the antenna, and the driving 
control torques are determined from the knowledge of the 
full antenna dynamics. Since only a small part of the an- 
tenna dynamics is measured (by encoders), the estimator 
is implemented to reconstruct the “missing” dynamics. 
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The controller and the estimator designs consist of ad- 
justing their gains through proper determination of the 
controller and estimator weights. For the antenna model 
of order n, 2 n 2 weights have to be determined (n 2 for the 
controller and n 2 for the estimator). This number is cus- 
tomarily dropped to 2n weights (n for the controller and 
n for the estimator). But in^spite of this drastic reduc- 
tion, the number of weights is still too large to make the 
search for the best weights reasonable (typically n = 40 
for the antenna, i.e., 80 weights have to be found). The 
difficulty arises because no general procedure for weight 
determination is available, and the known procedures deal 
with simplified and/or specific cases. 

The weight determination presented here becomes sim- 
ple due to several properties of the antenna and the com- 
pensator investigated in this article. First, it is shown that 
the tracking and flexible motion of the antenna are almost 
separated. The tracking part consists of four states (ele- 
vation and azimuth encoder readings and their integrals); 
thus, instead of dealing with a model of order n, one ob- 
tains two separate models of order 4 and n — 4. Secondly, 
the LQG compensator is balanced such that the controller 
and compensator efforts are the same. For the balanced 
compensator, the weights of the controller and estimator 
are the same, thus instead of 2n weights, n weights need to 
be determined. Thirdly, it is shown that each component 
of the balanced compensator (consisting of two states) is 
almost independent of others. Thus, weights for each com- 
ponent are determined separately. In consequence, the 
search for 2n weights becomes a series of searches for 2 
weights, which obviously is not a difficult task to perform. 

As mentioned, the PI controller is easy to implement 
due to its simplicity. But the implementation of the 40- 
state model-based LQG compensator is not an easy task. 
It would result in a complex algorithm, and would be a 
huge computational burden. Therefore, a simplification of 
the LQG compensator is an important implementation re- 
quirement. The simplification is obtained through order 
reduction of the compensator. The size of the compen- 
sator is reduced, but one has to find a reduction procedure 
such that the reduced-order compensator is stable and its 
performance is still close to the full-order compensator. 
This task is solved by introducing the pole mobility index. 
The pole mobility index characterizes the importance of 
the components of the balanced compensator. The states 
with small pole mobility index are truncated, and the trun- 
cation marginally affects the closed-loop dynamics. The 
closed-loop system with the reduced compensator is sta- 
ble, and its performance is close to the full-order compen- 
sator. It is shown that the reduced-order compensator of 
12 states is stable and its performance is close to the per- 
formance of the full-order compensator of 40 states. 


II. Problem Statement 

The closed-loop system with an LQG compensator is 
shown in Fig. 1, with the plant state-space triple ( A , B , C) y 
the process noise v of intensity V > and the measurement 
noise w of intensity W, where both v and w are uncorre- 
lated: 


V = E(vv t ) ) W = E(ww t ) 
E(vw t ) = 0, E(v) = 0, E(w) = 0 


( 1 ) 


where E( .) is an expectation operator. It is assumed that 
W — 1 without loss of generality. The task is to determine 
the controller gain ( K p ) and estimator gain (K e ) such that 
the performance index J 

( x T Qx + u T Ru)dt^ (2) 

is minimal, where R is a positive definite input weight ma- 
trix, and Q is a positive semidefinite state weight matrix. 
It is assumed that R — 7, also without loss of generality. 
The minimum of J is obtained for the feedback u = —K p x, 
where the gain matrix 

K p = B t S (3) 

is obtained from the solution S of the controller Riccati 
equation (CARE) [4,5]: 

A t S + SA- SBB t S + Q = 0 (4) 



The optimal estimator gain is given by 

K e = PC T (5) 

where P is the solution of the estimator Riccati equation 
(FARE): 


AP + PA t - PC T CP + V = 0 (6) 

The LQG compensator performance can be signifi- 
cantly modified through variations of weight Q and vari- 
ance V. Although V is formally predetermined by the 
process noise v , it can be modified in a search for a more 
suitable solution [5,6]. The determination of the weight 
and covariance is addressed in the following section. 
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Another nontrivial issue addressed here is the order of 
the compensator. Although the size of the plant deter- 
mines the size of the compensator, in many cases a full-size 
compensator is not acceptable for implementation due to 
its complexity. Thus, its order must be minimized in such 
a way that the reduced compensator maintains the stabil- 
ity and performance of the full-order compensator. The 
solution to this problem is found through the approximate 
balancing of the CARE/FARE equations. 


III. Quasi-Separation of Flexible and 
Tracking Subsystems 

The open-loop (or rate-loop) state-space representa- 
tion (A,B,C) of the DSS-13 antenna includes the input 
uj — [u pe u pa ], which consists of the rate commands 
in elevation ( u pe ) and azimuth (u pa ) and the outputs 
yj = [j/pe Vpe] and yj = [y ie yi a ], which consist of the 
elevation and azimuth angles ( y pe , y pa ) and their inte- 
grals (yie, yia)- Divide the state vector x of the open-loop 
antenna model into the tracking Xt and flexible xj parts 

X T = [xj xj] (7) 


where xj — [yj yj] and xj are the remaining states. It 
can be shown [2] that, in this case, the rate-loop represen- 
tation (A, B,C) has the form 


A = 


Ai A . tJ 1 , = 

'B t ' 

0 Af J ’ 

. B f . 


C=[Ct 0 ] 


( 8 ) 


and that 


II II < II 5 / II. \\A t f\\<\\A t \\, || A t/ 1| <|| A/ || (9) 

For the DSS-13 antenna, || B t || < 10“ 5 , || Bj || > 1, || 
A t j || < 10“ 3 , || Aj || > 10, and || A t || = 1. Thus, the 
states of the tracking part are much weaker than the states 
of the flexible part. The strong and weak signal flow is 
shown in the block diagram of Fig. 2. The strong states of 
the flexible subsystem and the weak states of the tracking 
subsystem are shown in Fig. 3, which presents the transfer 
function plots of the rate-loop systems due to elevation 
rate command. 

In the LQG design, the performance index is minimized 
and the minimum is obtained for K — B T S ) and S' is a 
solution of the Riccati Eq. (4). Divide S and K into parts 
related to the triple (A, B, C) in Eq. (8) 


5 = 



Sts' 
Sj J 


K = [K t K f ] 


( 10 ) 


so that Eq. (4) can be written as follows 


Aj S, + S t A t - S t B t Bj S t + Q t - A t , =0 (11a) 
AjS t j + S t jAj + S t A tJ - Kj K f = 0 (lib) 
Aj Sj + Sf Aj + SjBj Bj Sj-\-Qf— Aft — 0 (11c) 


where 


(12a) 


Ajt = AjjStj + SjfA t f - Sj f B t K f 
ASjBjBJ S t j 

A t ] = S t B t Bj Sj f + StfBjKt 
Kj = Bj S tJ + BjS f , K t = Bj S t + Bj Sjj (12b) 


Taking a closer look at Eqs. (12), notice that there exist 
weights Q t and Qj such that the gain I\j depends on the 
flexible subsystem only. Namely, for large enough Qj , such 
that || Qj || > || A f t ||, the solution Sj of Eq. (11c) is 
independent of the tracking system, and for small Q t) one 
obtains WBfStjW <C \\BjSj\\. In terms of Eq. (12a), the 
latter inequality means that the gain Kj depends only on 
the flexible subsystem. However, due to the master-slave 
relationship between flexible and tracking subsystems, the 
situation is not quite symmetric: There are no such Q t 
and Qj that the gain K t depends only on the tracking 
subsystem. To understand this, note that the term “small” 
has a different meaning for Qj and Q t . Magnitudes of 
small Qj and small Q t are of different order, namely small 
Qj is such that Qj < 10 -7 , and small Qt is such that 
Q t < 1. Therefore, by increasing Q t in order to obtain 
|| Qt || ^ || A. t / || one obtains || Bj Sj t || and || Bj S t || of the 
same magnitude. According to Eq. (12b), the latter fact 
means that the gain K t depends on the flexible subsystem, 
as well as on the tracking subsystem, and the solution St 
of Eq. (11a) is dependent on the flexible subsystem. This 
property can be validated by observation of the closed- 
loop transfer functions for different weights (Fig. 4). The 
plots show that the variations of Qj changed the properties 
of the flexible subsystem only, while the variations of Q t 
changed the properties of both subsystems. 

The weight Q t should be large enough to achieve the 
pointing performance requirements, and the increase of Q t 
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causes increasing dependency of the gains on the track- 
ing system. For this reason the above independence be- 
comes a quasi-independence in the final stage of controller 
design. Nevertheless, the separation in the initial stages 
of controller design is very strong. The design consists, 
therefore, of the initial choice of relatively small weights 
for the tracking subsystem and determination of the con- 
troller gains of the flexible subsystem. It is followed by 
adjustment of weights of the tracking subsystem and a fi- 
nal tuning of the flexible weights. 

The quasi-separation principle discussed above in the 
case of the controller design is also valid in the case of 
estimator design, since the compensator design consists of 
the independent designs of a controller and an estimator 
[4,5]. Additional properties of the LQG compensator that 
arise in controls of flexible structures are discussed in the 
following section. 


of Q = C T C and V = BB T . Gawronski [11] has shown 
that the balanced LQG representation exists in the case of 
general Q and V, and has derived the transformation to 
the balanced LQG representation. 

Let (A, B y C ) be a state-space triple of an open-loop 
system. Its controllability and observability grammians 
W c and W 0 are positive-definite and satisfy the Lyapunov 
equations 


AW C + W c A t + BB t = 0 
A t W 0 + W 0 A + CC T = 0 


(14) 


The system representation is balanced in the sense of 
Moore [12] if its controllability and observability gram- 
mians are diagonal and equal: 


IV. Balanced LQG Compensator for 
Flexible Structures 

In this section the flexible subsystem is considered only 
(subscript “/” is dropped in this section for simplicity of 
notation). A flexible structure is defined as a controllable 
and observable linear system with distinct complex con- 
jugate pairs of poles {N poles, N is even) and with small 
real parts of the poles. In other words, it is a linear sys- 
tem with vibrational properties. In the Moore balanced 
coordinates, it consists of n = Nf 2 components [7,8], and 
each component consists of two states. 

An approximately balanced LQG compensator is con- 
sidered. An approximate equality between two variables is 
used in the following sense: Two variables x and y are ap- 
proximately equal (x = y) if x — y-f e, and ||e||/ \\y\\ 1. 

The block diagram of a closed-loop flexible system with 
the LQG compensator is shown in Fig. 5. Similar to the 
balancing of controllability and observability grammians 
is the balancing of CARE and FARE equations. Namely, 
there exists a diagonal positive definite M = diag (//;), i = 
1, n, /i* > 0, such that 

S — P ~ M (13) 

A state-space representation with the condition 
Eq. (13) satisfied is called an LQG balanced representa- 
tion, and //*, i — 1, ...,n represents the characteristic val- 
ues of Jonckheere and Silverman [9] and Op- 

denacker and Jonckheere [10] have shown that a balanced 
solution for CARE and FARE equations exists in the case 


W C = W 0 = T 2 } T = diag(y u ...,y N ), i=l,...,tf (15) 

where ji > 0 is the ith Hankel singular value of the system. 
In [11] it is shown that, for flexible structures, the balanced 
representation (in the Moore sense) produces diagonally 
dominant solutions of CARE and FARE, and in the case 
of Q = V produces approximate LQG balanced solutions 
S and P, such that S =. P M. Assume a diagonal 
weight matrix Q : 

Q = diag(q i I 2 ) ) i=l,... f n (16) 

Then there exists ^ < q Q { , where >0, i = 1, n, such 
that S = diag (sil 2 ) is the solution of Eq. (4), where 

* = (fa - l)/27, ? . fit = 1 + 2*7?/C<w< (17) 

The proof is presented in the Appendix. The plots of S{ 
with respect to q t and 7 ,• are shown in Fig. 6. They show an 
increase of s, with the increase of weight q t , and decrease 
of Si with increase of 7? or Ct^i* 

A similar result is obtained for the FARE equation, 
namely, for a diagonal V : 

V = diag(vili) t i— l,...,n (18) 

there exists Vi < v 0 { , where v Q i >0, i — 1, such that 
P = diag(pil 2 ) is the solution of Eq. (6), where 
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Pi = (A< - l)/27, ? , = 1 + 2«*7?/<<« (19) 


Ppi — Cci / C* — Ctoil&ci 


(24) 


If the ith diagonal entry of P and the respective entry 
of S are equal, say to //*, i.e., 

Vi = Si - Vi ( 20 ) 

the zth component is LQG balanced; if S and P are equal 
the system is LQG balanced. If 5, P, and M are diago- 
nally dominant, i.e. v* + = s, + e 3l - — Vi> with e v * 

and small (| t vi /vi | < 1, | | <C 1), the system is 

approximately LQG balanced. 

It follows from Eqs. (17) and (19) that for Q = 
diag(qi) = V = diagfa), the system is approximately LQG 
balanced, such that 

S — P — M = diag(vi) 

( 21 ) 

/X, = (A - l)/27?, A 2 = 1 + 2 9iT 2 /C> w 


i.e., is a ratio of closed- and open-loop damping factors, 
or it is a ratio of open- and closed-loop resonant peaks. 
Therefore, if a suppression of the ith resonant peak by the 
factor P p i is required, the appropriate weight ^ is deter- 
mined from Eq. (17): 

q, = 0.5(/?p,- - l)CiW<7i" a ( 25 ) 

An alternative interpretation of Pi is as a ratio of the 
open-loop Hankel singular value to the closed-loop Hankel 
singular value: 

A = ilihli ( 26a ) 

or a ratio of variances of open-loop (<x^) and closed-loop 
(<r^) states excited by the white-noise input [13]: 

A = (26b) 


Next it is shown that the weight Q 

Q - diag{ 0, 0, ... ,qih , ---0, 0), qi < q 0 i (22) 

shifts the ith pair of complex poles of flexible structure, 
and leaves the remaining pairs of poles almost unchanged. 
Only the real part of the pair of poles is changed (just mov- 
ing the pole apart from the imaginary axis and stabilizing 
the system), and the imaginary part of the poles remains 
unchanged. Namely, for the weight Q as in Eq. (22), the 
closed-loop pair of flexible poles (A cr », ±j\ c u) relates to 
the open-loop poles (A or {, ±jA 0 ,i) as follows: 

(A cr j, dbjAcii) — (Ppi^orii ijA^i), i — (23) 

where P pi is defined in Eq. (17). The proof is presented in 
the Appendix. 

The real part of the poles is shifted by /? p », while the 
imaginary part remains unchanged. The above proposition 
has additional interpretations. Denote the real part of 
the open-loop pole by A ort - = — Ci w * an d the real part of 
the closed-loop pole by A cr ; = — Cc*u>*j n °t e also that the 
height of the open-loop resonant peak is a 0 i = 
where k is a constant, and the closed-loop resonant peak 
is Gt c i — /c/2£ C jU>,-. From Eq. (23), Ppi — A cr j/A or ,’; thus, 
one obtains 


The proof is presented in the Appendix. 

The plots of Ppi with respect to qi and 7, are shown in 
Fig. 7. They show relatively large P p i even for small qi , i.e., 
a significant pole shift to the left. Also, P p i increases with 
the increase of 7*, and decreases with the increase of 
i.e., there is a significant pole shift for highly observable 
and controllable states with small damping. In terms of 
the transfer function profile, the weight qi suppresses the 
resonant peak at frequency w* while leaving the natural 
frequency unchanged. Due to weak coupling between the 
states, the assignment of one pair of states does not signif- 
icantly impact other states. Thus, the weight assignment 
can be done separately for each pair of states. 

The estimator poles are shifted in a similar manner. 
Denote 

V = diafli(0,0, ...,Ui/2,...0,0), < v oi (27) 

then for the weight V as in Eq. (27) and v t - < v 0 i y the 
estimator pair of poles (A er i, hj\ e a) relates to the open- 
loop poles (A or {, ±jA 0 u) as follows: 

(A er j, ±j*A e ii) = {Pei^orit ijA 0 ii), i— (28) 

where /? ei * is defined in Eq. (19). 
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(29) 


The limiting values q 0 i and v 0 i are determined. Their 
values are rather fuzzy numbers. Despite their fuzziness, 
they are not difficult to determine. There are several indi- 
cators that qi is approaching q 0 i , or that Vi is approaching 
v ot . For the controller, q 0 i is the weight for which the 
ith pair of complex poles of the plant departs from the 
horizontal trajectory in the root-locus plane, or it is the 
weight for which the ith resonant peak of the plant trans- 
fer function disappears (the peak is flattened). For the 
estimator, v 0 i is the covariance for which the ith pair of 
complex poles of the estimator departs from the horizon- 
tal trajectory in the root-locus plane, or it is a covariance 
for which the ith resonant peak of the estimator transfer 
function disappears. 

It is crucial from an implementation point of view to ob- 
tain a compensator of the smallest possible dimension that 
preserves the stability and performance of the full-order 
compensator. Although the size of a plant determines the 
size of a compensator, the plant model cannot be reduced 
excessively in advance in order to assure the quality of the 
closed-loop system. Therefore, compensator reduction is a 
part of the compensator design. The balanced LQG design 
procedure provides this opportunity. 

To successfully perform the compensator reduction, an 
index of the importance of each compensator component 
is introduced. In the open-loop case, Hankel singular val- 
ues serve as reduction indices. In the closed-loop case, the 
characteristic values of the system seem to be good can- 
didates for the reduction indices, as suggested in [9]. This 
is not a good choice, however, since the characteristic val- 
ues do not properly reflect the effectiveness of the com- 
pensator. The effectiveness of the closed-loop system can 
be evaluated by the relative suppression of the closed-loop 
output when compared to the open-loop output. Thus, the 
ratio of variances of the closed-loop output and the open- 
loop output excited by the white noise input is an appro- 
priate measure of the suppression (alias the compensator 
performance). It will be shown later that the suppression 
depends on pole mobility in the complex plane. Therefore, 
if a particular pair of poles is easily moved (i.e., in the case 
when small weight is required), the respective states are 
easy to control and to estimate. On the contrary, if a par- 
ticular pair of poles is difficult to move (i.e., a large weight 
is required to move the poles), the respective states are 
difficult to control and to estimate. In the latter case, the 
action of the compensator is irrelevant, and the states that 
are difficult to control and estimate can be reduced; thus, 
pole mobility is a good indicator of the importance of a 
particular compensator state. 

Consider an LQG balanced system, and denote the pole 
mobility index i r, by 


T,- - 0.5 (Pi - 1) 

Note that for # = 1, the s’th pole is stationary and i r t - is 
equal to zero; for a shifted pole, one obtains 0i > 1 and 
7 Ti > 0. If for small x,- a small pole shift (in plant, as well as 
in estimator) is observed, this component can be reduced. 

Another useful interpretation follows from Eqs. (29) 
and (26): 

* = 0.5( 7o 2 , - iDhli = 0.5 {o*t - (30) 

i.e., the pole mobility index is proportional to the relative 
change in the white noise response of the open- and closed- 
loop systems. Furthermore, from Eqs. (21) and (29), it can 
be shown that the pole mobility index is a product of the 
square of a Hankel singular value and the characteristic 
value of the system: 

*i = 7 1 m ( 31 ) 

Thus, 7 combines the system observability and controlla- 
bility properties of the open-loop system with the closed- 
loop performance. The more heavily weighted the com- 
ponent, the larger its pole mobility index, see Fig. 8(a). 
Also, the larger the Hankel singular value of the compo- 
nent, the larger the corresponding pole mobility index, c.f. 

Fig. 8(b). 

The matrix n of pole mobility indices is defined as 

II = c?iafli(7r 1 ,ir2,...,7r n _j,x n ) (32) 

and it is obtained from Eqs. (15) and (29) as 

n = T 2 M (33) 

In the following, a reduction technique is discussed. As- 
sume n in Eq. (32) has a descending order, i.e., x,- > 0, 
7T t * + i < x*, i = and divide n as follows: 

H = diag(U r} n*) (34) 

where n r consists of the first k entries of n, and n* the 
remaining ones. If the entries of U t are small in compari- 
son with the entries of n r , the compensator is reduced by 
truncating its last n — k states. 
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(35a) 


It is shown in [11] that a system with the reduced-order 
compensator obtained by reducing states with small II* is 
expected to be stable. That is, although it is not guar- 
anteed, there is a well-founded expectation to obtain a 
stable system with the reduced-order compensator. Also, 
the estimation errors of the full-order and reduced-order 
compensators are approximately the same. 


V. Closed-Loop System 

The LQG compensator configuration for the DSN an- 
tenna is shown in Fig. 1. The tracking command y c is 
compared with the estimated antenna position y p , and the 
error e — y p — y c along with the integral e* of the error 
e are the plant inputs. The equations for the integrator, 
plant, and the estimator are, respectively, 


ki — e 

x = Ax + Bu, y p = C p x , xj = Cjx 
x = Ax + Bu + K e (y p — C p x) 
u ~ —KjXf — K p e — KiC{ 

Denoting xj f = [ej x T x T ], one obtains 

Xd = Acix c i + B c iy c , y = C c \x c \ 

where 


(35b) 


(35c) 


(36a) 



■ 0 

0 

c p 


- -i - 



A d = 

-BKi 

A 

—BKjCj - BK p Cp 

, B c i = 

BK p 

, Cel = [0 Cp 0] 

(36b) 



K e C p 

A - BK/Cj - BI<pC p - K e Cp . 


-BK P . 




The triple (A c j, 1 B c i, Cel) for the LQG system with the reduced-order estimator (A r , B r) C r ) is as follows: 



o 

o 

■p 

r 


■ -i - 

A c i — 

-BI<i A -BKjCjr - BKpCpr 

^ B c i — 

BKp 


_ — B r I\i I\. er C p A. r — B r I\ fCfr — B T K p Cpr d\ er C p r - 


-B r Kp. 


, C c i = [0 C p 0] (37) 


where C pr and Cj r are obtained from the partition of CV, C r = [C pr Cj r ]• 


VI. Compensator Design Algorithm 

Weights for the balanced LQG controller and estimator 
are identical. Therefore, in the algorithm, the controller 
and estimator gains are adjusted simultaneously. The pro- 
cedure for the antenna LQG compensator design is a se- 
quential one. First, for the ad hoc (but relatively small) 
chosen tracking subsystem weights, the flexible subsystem 
weights are determined (recall that the controller and es- 
timator weights are the same). Second, the adjustment 
of the tracking system weights is performed, followed by 
final adjustment of the flexible system weights. The flexi- 
ble subsystem weights are determined sequentially, which 
gives more insight into the system performance and sim- 
plifies the procedure. 


The estimator order is determined as a part of the 
weight tuning process. Only the modes with large index 
Xj are considered. If the number of flexible modes is ny, 
the number of disregarded modes is n 0 , and the size of 
the tracking system is n t , then the controller order n c is 
n c — n t + 2 (nj — n 0 ). The following LQG compensator 
design algorithm is proposed: 

(1) Determine the plant state-space representation in 
the form of Eq. (8), consisting of flexible and track- 
ing parts. 

(2) Choose ad hoc but reasonably small weights and 
variances for the tracking part Q t — V t = Q tah- 
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(3) For each balanced coordinate of the flexible part, 

choose the weight ft, i = and define the 

weight and covariance matrices Q/* = 

diag (0,0,..., ft, ft,0,0,...0), V fi = <?/* so that 
the closed-loop system performance for the weight 
Qi = Q/») and the covariance V* = Q; is 

optimized. For example, determine the weights ft to 
impose the required pole shift or to suppress the *th 
resonant peak to the required level without depre- 
ciating other properties of the closed-loop transfer 
function. Note that for small ft, only the it h pair 
of poles is shifted (to the left), and the remaining 
poles are almost unaffected. Disregard the modes 
with small index i r, for which the weighting does 
not improve the closed-loop system performance. 

(4) For the already determined Qj and Vj , tune weight 
Qt and assume covariance V t = Qt to obtain im- 
provements in tracking properties of the antenna. 

(5) Adjust flexible subsystem weights, if necessary. 

VII. Applications 

A balanced LQG compensator was designed for the 
DSS-13 antenna. The DSS-13 antenna model consists of 
two tracking states (azimuth and elevation angle) and 13 
flexible modes (or 26 balanced states). The preliminary 
weights qu = q pe = ft a = q P a = 1 for the tracking sub- 
system (for ft and y p ) and zero weights for the flexible 
subsystem (ft = ft = ••• — $13 — 0 f° r a H 13 modes) 
were chosen. The closed-loop system step response is pre- 
sented in Fig. 9 (azimuth encoder reading due to azimuth 
command) and the magnitudes of the closed-loop trans- 
fer function in Fig. 10. Both figures show that flexible 
motion of the antenna is excessive and should be damped 
out. This is achieved by adjusting weights for the flexible 
subsystem. For the tracking weights as before, the weight 
for the first mode (2.32 Hz) is chosen to be ft — 10" 7 , and 
the remaining weights are zero, obtaining the closed-loop 
system responses as shown in Figs. 11 and 12. One can see 
that the 2.32-Hz resonance peak in the azimuth command 
response (Fig. 11) has disappeared, along with most of the 
flexible motion in the azimuth step response (Fig. 12). The 
elevation motion is unaffected, however, since the azimuth 
mode is almost nonexistent in the elevation motion. 

The weight should be chosen carefully. Weight that is 
too small (e.g., 3 x 10“ 9 in the case considered) will not 
suppress the resonant peak, Fig. 13(a). Weight that is 
too large (e.g., 10" 5 ) will deteriorate the tracking perfor- 
mance: for the overweighted mode the transfer function is 
pressed down within a wide frequency range, Fig. 13(b). 
The proper weight suppresses the resonant peak, leaving 
the other peaks unchanged. 


Similar procedures have been applied for the second 
(2.64-Hz), third (4.26-Hz), fourth (3.77-Hz), fifth (7.92- 
Hz), sixth (4.47-IIz), seventh (3.38- Hz), eighth (5.98-Hz), 
ninth (7.32-Hz), and tenth (9.48-Hz) modes, with weight 
10“ 7 for each mode. As a result, the suppression of the 
remaining flexible motion is observed as shown in Fig. 14. 
Weights for the remaining modes (eleventh through thir- 
teenth) have been set to zero. 

The root locus of the closed-loop system due to weight 
variations of the 7.92-Hz (fifth) mode is shown in Fig. 15. 
The figure shows the horizontal departure of poles to the 
left (stabilizing property). It confirms the properties of 
the weighted LQG design described previously. 

In the next step, the tracking properties of the system 
are improved by proper weight setting of the tracking sub- 
system. Namely, setting the integral weight to ft e = q ia = 
70 and the proportional weight to q pe = q pa = 100 im- 
proves the system tracking properties, as shown in Fig. 16 
(small overshoot and settling time) and in Fig. 17 (ex- 
tended bandwidth up to 2 Hz). However, by improving the 
tracking properties, the transfer function has been raised 
dramatically in the frequency region of 1 to 3 Hz, which 
forces the first two modes located in this region to ap- 
pear again in the step response. By sacrificing a bit of 
the tracking properties, the flexible motion in the step re- 
sponse is removed. This is done by slightly increasing the 
weights of the flexible subsystem, setting them as follows: 
?i = ?2 = 73 = 94 = 95 = 96 = 10 -6 , <77 = Qs = 10~ 7 , and 
q 9 = fto — 10" 5 . The closed-loop system response with 
satisfactory tracking performance is shown in Figs. 18 and 
19 (small overshoot, settling time, and 1-Hz bandwidth 
are observed). 

The reduced-order compensator is obtained through 
evaluation of pole mobility indices i r*. The plot of is 
shown in Fig. 20. Reducing the order of the estimator 
to 12 states [the first four are tracking states (states 1 to 
4), the next six are flexible states (states 5 to 10), and 
the last two states are nonflexible components of the rate 
loop model (states 25 and 27)] gives the stable and accu- 
rate closed- loop system. The reduced-order compensator 
is compared with the full-order compensator in the step 
response plots in Fig. 21, and the transfer function plots 
in Fig. 22, showing satisfactory approximation. 

VIII. Conclusions 

A method for designing a reduced-order compensator 
for the DSS-13 antenna has been presented. - A balanced 
LQG compensator has been introduced that uses the same 
amount of effort to control and to estimate the system. 
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The properties of the balanced LQG system are used 
to obtain a reduced-order compensator for the antenna. 
This compensator preserves the stability and performance 


of the full-order compensator. The performance of the 
reduced-order compensator has been verified by simula- 
tions. 
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Fig. 3. Magnitudes of transfer function of tracking and flexible subsystems for the eleva- 
tion rate input. 
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Fig. 4. Magnitudes of transfer function of the closed loop system for different weights of 
elevation angle to elevation command. 


93 


i 1 1 1 1 ifl 





94 


















HANKEL SINGULAR VALUE y 


Fig. 7. Pole shift factor /3 versus: (a) weight q and (b) Hankel singular value *y. 







HANKEL SINGULAR VALUE y 


Fig. 8. Pole mobility Index tt versus: (a) weight q and (b) Hankel singular value 7, 
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I. 9. Azimuth encoder response to azimuth step command for the unit proportional and in- 
tegral weights In azimuth and elevation, and zero weight for flexible subsystem. 



FREQUENCY, Hz 
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Fig. 10. Transfer function (elevation and azimuth encoder due to azimuth command) for 
weights the same as those In Fig. 9. 
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Fig. 11. Transfer function (elevation and azimuth encoder to azimuth command) for weights 
the same as those in Fig. 9, but for qn = 10“ 7 . 
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Fig. 12. Azimuth encoder response to azimuth step command for weights the same 

as those in Fig. 11. 
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Fig. 13. The first component: (a) underweighted and (b) overweighted. 
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TIME, sec 

Fig. 14. Azimuth encoder response to azimuth step command for weights the same as 
those In Fig. 9, but for qn = ... = <7/10 = 10 -7 , 



REAL PART 

Fig. 15. Root locus for 7.92-Hz mode versus weight. 




Fig. 16. Closed-loop response to step Input for proportional weight 100 and Integral weight 
70 (both in azimuth and elevation), and weights for flexible subsystem equal to 10 -7 (azimuth 
encoder to azimuth command). 



Fig. 17. Closed-loop transfer function for weights the same as those In Fig. 16 (elevation en- 
coder and azimuth encoder to azimuth command). 
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1.2 



TIME, sec 

Fig. 18. Azimuth encoder response to azimuth step command for proportional weight 100, 
integral weight 70 (In azimuth and elevation), and for flexible weights qn ■ = g/6 = 10 6 , 

q n = q f a = 10~ 7 , and q f g = gno = ™“ 5 - 
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Fig. 19. Closed-loop transfer functions for weights the same as those In Fig. 18 (elevation 
encoder and azimuth encoder to azimuth command). 
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Fig. 22. Transfer function (azimuth encoder to azimuth command) of the DSS-13 antenna with 
full- and reduced-order compensators. 



Appendix 


I. Selected Properties of Flexible Systems 

The balanced grammian for a flexible system with n 
components (or N — 2n states) has the following form 
[7,8]: 

r S diagiji , 71 . 72 , 72 , 7„ , y n ) (A-l) 

The system matrix A is almost block diagonal, with dom- 
inant 2x2 blocks on the main diagonal: 


There are two solutions of Eq. (A-6), but for a stable sys- 
tem and for g* = 0 it is required that s, = 0, therefore 
Eq. (17) is the unique solution of Eq. (A-6). 

III. Proof of Eq. (23) 

For small qi , the matrix A of the closed-loop system is 
diagonally dominant A 0 — diag (A 0 ;), i = 1 , n, and 
A 0 i = Ai - BiBj Si. Introducing Eq. (A-3), one obtains 

A 0 i — A{ + 2s,- 7^ (A{ + Aj) (A-7) 


A 2* diag(Ai),Ai 


-Qu>i -u>i 

Wj — CiW* 


i= (A-2) 


and introducing Ai as in Eq. (A-2) to Eq. (A-7) one obtains 


where w,- is the zth natural frequency of the structure, and 
is the zth modal damping. The matrices B and C are 
divided into two blocks, comparably to A } B T = [Bf , 
tfj,..., and C = [Ci,C 2 , C n ], with the following 

property: 


A 0 i — 


-PpiCiVi 

-/?p*C«w* 


-'ptS 

U>i 


with f3 p i as in Eq. (17). 


(A-8) 


BiBj Sf CjCi “ + Af ) (A-3) 

II. Proof of Eq. (17) 

Due to the diagonally dominant matrix A for a flexi- 
ble structure in balanced representation, and for Q as in 
Eq. (16), there exists qi < q 0 i, i = 1, n } such that the so- 
lution S of the Riccati Eq. (4) is also diagonally dominant 
with 2x2 blocks Si on the main diagonal: 

Si “ Sil 2 , Si > 0, » = 1, n (A-4) 

Thus, Eq. (4) turns into a set of the following equations: 

s,(Ai + Aj) - sfBiBj + qj 2 = 0, i = 1, .... n (A-5) 

For a balanced system BiBj — —yf (Ai + Aj) and A, + 
Aj = -2CiUih, see Eqs. (A-3) and (A-2), respectively. 
Therefore, Eq. (A-5) is now 

s 2 + s > hi ~ O.bqi/Ci^iJi — 0. * = 1. •••> n (A-6) 


IV. Proof of Eqs. (26a) and (26b) 

In order to prove Eqs. (26a) and (26b), the closed-loop 
Lyapunov equation is considered: 

(A - BB t S)T 2 + r^(A - BB T S) T + BB T ~ 0 (A-9a) 

or, for the zth pair of variables, 

(Ai-BiBjsi)y 2 ci +y 2 cl (A t -BiBj Si ) T +BiBj 9£ 0 (A-9b) 

Introducing Eq. (A-3) gives 

y 2 ci + 2yW 0iS i-y 2 0i 2 0 (A-10) 

or 

llihl S£ 1 + 28^ (A-ll) 

Comparing Eq. (A-ll) and Eq. (17) gives Eq. (26a). 
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